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ABSTRACT 

We present a cascading model for a single jet, inspired to QCD 
and to the phenomenological analysis of multiplicity distributions. 
The model, describing as it does a two dimensional evolution in 
virtuality and rapidity, allows analytical predictions for clan analysis 
parameters to be made. 



Talk presented by R. Ugoccioni 
at the XXV International Symposium 
on Multiparticle Dynamics, 
Stara Lesna, Slovakia, 12-16 Sept. 1995. 
To be published in the Proceedings. 



1. Introduction 
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Since its proposala, clan analysis has been widely used experimentally in order 
to interpret the negative binomial (NB) regularity and as a tool to classify different 
reactionsi. It is especially important to stress the presence of the NB regularity 
when the multiplicity distributions (MD) are studied in intervals of rapidity, where 
conservation laws don't play a dominant role and therefore an understanding of the 
dynamics is most interesting and desirable. 

Recently clan analysis has been extended in three directions: the regularity is 
better satisfied in the case of single jets, as is suggested by Monte Carlo studiesi 
and experimentally by the observation that separating the 2-, 3- and 4-jet sample 
of events restores the regularity that was violated (in the shape of the multiplicit 
distribution, but not in the clan structure parameters) in the full sample of event 
Secondly, one should remember that, in going from the hadronic level to the partonic 
one via generalized local parton-hadron dualityi, the NB regularity is better satisfied. 
This suggests that an explanation of the regularity should be sought at partonic level. 
The third extension comes by recognizing that the clan interpretation of the negative 
binomial distribution (NBD) is actually more general than the NBD itself, and leads 
to the idea that the emission of partons in a single jet is a two-step process!. 

The difficulty of analytical calculations in QCD, and a guess on the real nature of 
clans in multiparticle production (they might very well be true physical objects) has 
prompted us to develop an analytical parton shower model for a single jet inspired to 
QCD and at the same time close to the phenomenological observations. 

In this paper we first discuss some properties of two-step processes in general, and 
once this framework has been established, we describe the model and its results. 



2. Two-step processes 

Assume that the parton production process involves two independent steps: in the 
first step N objects {N = 0, 1 . . .), (which we call clans, and the clans known from the 
NBD are now a particular case), are produced with probability pn and generating 
function f{z) (in all quantities in this Section a dependence on the jet energy is 
understood). In the second step, each clan produces partons with probability g„. 
(with i labeling the clans: i = 1,2 . . . N) and generating function g{z). Since clans are 
identified by the final partons they generate, requiring that each clan produce at least 
one parton (i.e, go = 0) makes the number of clans unambiguously defined. Notice 
that we do not assume that all clans are created equal: the multiplicity distribution 
(MD) of an individual clan can depend on a set of parameters as qmiO'- if the value 
of ^ for a clan is independent of that of other clans and on the number of clans, one 
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can define an average clan: 

qn.= jqnXmm. (1) 

where 0(^) is the p.d.f. for producing a clan with that values of the parameters. 
Indeed this property will be used in Sec. 4. 

With these assumptions, the final partons MD, has generating function 

oo 

F{z) ^Y.^''Pn = fi 9{z) ) . (2) 

n=0 

If we turn now to look at the production in intervals of phase space, it is clear that the 
only dependence on the particular interval belongs to the second step, which is the 
step that produces the final partons. Choosing for definiteness to look at a rapidity 
interval which will be denoted Ay, we find that each clan produces partons according 
to 

oo 

g{z, Ay) = z^qniAy) , qoiAy) ^ (3) 

n=0 

where the fact that it is possible that a clan produces zero partons inside Ay has 
been emphasized (see also Fig. |1]). Obviously the final distribution is then given by 

F{z,Ay) = f{g{z,Ay)). (4) 

The parameters of interest in this paper are then the average number of clans, iV, 
and the average number of partons per clan, n^f^: 



-(o)^A..A - dgiz,^y) 



, dz 

z=l 



z=l 



n{Ay) 



dz 

where n{Ay) is the average number of final partons from Eq. @. 

2.1. Binomial convolution 

f{z) and g{z, Ay) can be redefined so that only clans which produce at least one 
parton in the interval Ay are considered: one simply subtracts qo{Ay) and rescales 
the distribution correspondingly: 

/ A ^ Ay) - qojAy) 

gi[z, Ay) = — — . 6 

1 - go (Ay) 

This implies then 

h{z) = f{[l-qo{Ay)]z + qo{Ay)). (7) 

It can be shown that this is equivalent to the condition that the probability that m 
of the N produced clans generate at least one parton in Ay is given by a binomial 
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distribution of parameter [1 — qQ{Ay)]. In particular, the average number of clans 
contributing to the interval Ay is in general 

N{Ay) = [1 - qo{Ay)]N (8) 

and therefore one obtains 

^""^^ NiAy) [l-go(Ay)]- 
Example lA. Suppose that f{z) is a Poisson distribution: 

f{z) = e^(^-^); (10) 

then the redefined distribution is still a Poisson distribution with parameter given by 
Eq. (@): 

F{z,Ay)=e^^{N{Ay)[g,{z,Ay)-l]}. (11) 
Example IB. Consider the case in which f{z) is a shifted Poisson distribution: 

f{z) = ^e^(^-i); (12) 

(such will be encountered in Sec. 3). Then the result is more complex: 

F{z,Ay) = {[1 -go(A|/)]^i(^, A|/) + go(A|/)}e^(^^)[^^(^'^^)-i] . (13) 

2.2. Compound Poisson distribution 

Alternatively, one can try and redefine the generating functions so that the final 
MD is a compound Poisson distribution; this is an implicit assumption made when 
one makes a fit with a NBD. We want to solve: 

F{z,Ay) = exp{A(A?/)[^2(^,A?/)-l]} (14) 
(?2(0,A|/) = (15) 

for \{Ay) and g2{z,Ay) and we obtain: 

X{Ay) = -log[F(0,Ay)] (16) 

^ , log[F{z,Ay)] 
92iz,Ay) = 1 + ^^^^^ . (17) 

Here g2{z, Ay) is a true probability generating function {d"'g2/dz"' > Vn) if and only 
if all the combinants of the distribution F are positive. In other words, in some cases 
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it may not be possible to carry out this redefinition. In case it is, the parameters of 
interest are A (Ay) as derived above and 



_ n(Ay) Nnf>( Ay) 

Example 2A. For a Poissonian /(z) we find the same result as in example lA: 

A(Ay) = [l-go(Ay)]iV (19) 
( A ^ ^(z,Ay) -go(Ay) 
*<-^^^> = l-*(Ay) l^"* 

SO that the two transformations coincide. In particular, A(A?/) is equal to the average 
number of clans in the interval Ay. 

Example 2B. For a shifted Poisson, on the other hand, they do not coincide: 

\{Ay) = [l-go(Ay)]iV-log[go(A2/)] (21) 

g2{z,Ay) = 1+ ' ^' J. . ^^^^ ^'^ 22 

A(Ay) 

and A (Ay) is not equal to N{Ay). In particular it should be noted that in full phase 
space (fps), iV(fps) = N but A(fps) goes to infinity. 



3. The GSPS model 



The GSPS model is introduced in order to solve analytically a two dimensional 
parton evolution in a single jeti'i. It describes an ancestor parton, which degrades in 
virtuality, and which we follow in rapidity, emitting clans of different virtuality and 
rapidity; each clan emits partons in a two dimensional cascade process, see Fig. |1]. 

3.1. Step One 

In Fig. |l| the first step is shown with thick lines. Because clans are by assumption 
emitted independently in a cascade fashion, we take the two branches in each splitting 
of the ancestor to be independently regulated by the same splitting function, for which 
a form inspired by QCD, suitably normalized by a Sudakov form factor term, has been 
chosen: 

PA{Qo\Q)dQo = pM)CA{Q)dQo = Ca{Q)^^ (cuk^) "^^^ ^ 

(23) 

The parameter A > thus controls the branching of the ancestor: to a larger value 
of A corresponds more branching. 



logQo 
logg 
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clan 1 



clan 2 clan 3 



Dy 

Fig. 1. Schematic representation of the GSPS model. Thick hnes indicate clans creation (step 1) 
and thin lines indicate cascading into final partons (step 2). Notice that only clan 1 and 2 produce 
partons in the interval Dy shown. 



In order to apply this factorization, we have to allow for local violations of the 
energy-momentum conservation law, still requiring its global validity, i.e., offspring 
partons of virtualities Qi can fluctuate according to: 



Qo + Qi^Q , IGeY <Qi<Q [i=0,l] . 



(24) 



Of course, constraints on the energy fraction carried away by daughter partons are 
also no longer valid, i.e., Zq + Zi ^ 1, and kinematic bounds in rapidity are loosened: 



|yi-y|<log^ [i=0,l] 



(25) 



This new condition modifies also the splitting kernel in z which now is decoupled: 

dz(j dz\ 



P{zo, Zijdzodzi oc 



Zq Z\ 



(26) 



Notice that, for each branch, we have taken the singular part of the AltareUi-Parisi 
kernel. 

Finally, let us mention that the rapidity of the first splitting is obtained from the 
energy W and the virtuality Q by the simple kinematic relationship: 



y = tanh-^ ^Jl-Q'^/W'^. 



(27) 



3.2. Step Two 

At this point, clans have been emitted with definite initial virtualities and ra- 
pidities. The simplest and most economic assumption is now that inside each clan 
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a cascade process develops in a way very similar to what has been described for the 
first step: each parton branches into two independent partons. The same form of 
splitting functions will be applied as in Sec. 3.1, but with a different parameter: 

Pa{Qo\Q)dQo = d (j^^ ; (28) 

to a larger value of a corresponds then a larger number of branchings, and therefore 
of final partons. The rapidity part is treated according to Eq. ( P3| ) and i^^- 

Finally, one should note that, if energy conservation were strictly enforced, a 
parton with virtuality less than 2 GeV could not split. Because of Eq. (^^ this is no 
longer necessary, but we keep it nonetheless as our cut-off procedure. 



4. The structure of the calculation 

The calculation has been described in detail in 0, here we will only outline its 
structure. 

Since, according to the model, clans of given virtuality and rapidity are produced, 
one must first calculate the generating function for the MD of a clan with initial 
virtuality Q and rapidity y to produce partons in the interval Ay, which we call 
g{z, Ay,Q,y). It would correspond to Eq. (^ in Sect. 2. It can be done based on 
step two of the model only. 

We then calculate the probability to emit a clan with initial virtuality Q and 
initial rapidity y, which we call (f){Q,y). This can be done based on step one of the 
model only. 

Finally one defines the generating function for the MD of an average clan similarly 
to what is done in Eq. ([l|): 

rW r 

g{z, Ay) = dQj dy g{z, Ay, Q, y)4>{Q, y). (29) 

This distribution is then to be convoluted with the generating function for the number 
of clans, f{z, W), which depends of course only on step one of the model. 

Example. It is perhaps best to illustrate this procedure by an example, which for 
simplicity only regards full phase spacei. Let us examine a pure birth model for a 
single clan: the MD is a shifted geometrical distribution 

g{z, fps, Q, y) = — — (30) 

z - i^{Q){z - 1) 

where the average value i^(Q) is assumed (for the sake of this example only: in the 
GSPS model its explicit form can be calculated analytically) to be a specific function 
of the clan virtuality Q: 

(P{Q)dQ oc — . (31) 

V 
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The MD of an average clan is given by 

^ J_ z du' ^ log(l - bz) 

log uJi z- u'{z - 1) u' log(l - 6) ^ ' 

where h = 1 — 1/v. It is a logarithmic distribution which, when convoluted with a 
Poisson distribution for the first step, gives a NBD. It should be pointed out that 
this example does not apply directly to the GSPS: as shown in the next section, the 
MD of an average clan is not a logarithmic distribution (but resembles it); the MD 
for clans is not a Poisson distribution, but a shifted Poisson distribution. 

5. Results and discussion 

The model can be solved analytically under a few mild mathematical approxi- 
mations: the full solution in rapidity intervals is however far from being simple and 
compact. The results are therefore presented in graphical form, but in order to get 
the general flavour it is interesting to quote the analytical result in full phase space: 
the MD of final clans is given by 

f{z, W) = ze^^[[N{W) - - 1]} ; (33) 

because at least one parton (the ancestor) is always present in the cascade, the total 
distribution of clans turns out to be a shifted Poisson distribution. The distribution 
inside an average clan is given by 

g{z, W) = [z + ^Ua{W) - ^ log [z+il- z) exp{ua{W)}]^ . (34) 



where 



rW (\o0:W\ 
Ua{W)^j^ p^{Q\Q)dQ = a\og{-^\ . (35) 

The final partons distribution is then obtained by inserting Eq. ( p4|) into Eq. (|33|) 
according to Eq. @). 

The clan parameters can then be calculated to be 

N{W) = l + A\og{^^\ (36) 



and 

(\ogW' 



nc{W) 



N{W) 




log 2 



(37) 



As expected, the average number of clans depends only on the parameter A which 
regulates the first step of the emission process; the average number of partons in an 
average clan depends on both parameters (whereas the distribution in a single clan 



7 



2 4 6 0.2 0.4 0.6 0.8 1 



Fig. 2. Left: Average number of partons in the shower, h{yc,W), as a function of the width of 
the rapidity interval yc obtained analytically in the GSPS model with A = 2, a = 1 at different 
maximum allowed virtualities = 50 GeV (solid line), 100 GeV (dashed line) and 500 GeV (dotted 
line). Right: The same quantity is plotted normalized to the value in full phase space as a function 
of the rescaled rapidity interval. 

of definite virtuality depends of course only on a). It should be noticed that in the 
limit of high initial virtuality (energy), W oo, the dependence of ficiW) on A 
disappears. 

The solution in rapidity intervals is, as explained in the previous section, a con- 
volution of Eq. (^) with the MD of an average clan; however, since we are here 
interested only on clan parameters properties, we will calculate only these, and not 
the full distribution. 

The actual calculations in rapidity intervals have been carried out for A = 2 
and a = 1, values which avoid nasty inessential calculations and make possible the 
analytical solution of the model without hurting its logic. 

In Fig. we plot the average number of final partons in both the standard form 
(on the left) and in a rescaled form (to the right), showing an interesting scaling 
behaviour with energy. This scaling in W is found to depend on the parameter a, as 
different values of a give different scaling curves, differently from the scaling found ini 
for N{Ay, W) / iV(fps, W) which by its own definition is independent of the mechanism 
at work inside clans. 

In Fig. 1^ we plot the parameter X{Ay,W), defined by requiring that full distri- 
bution be a compound Poisson distribution, and the corresponding average number 
of partons per clan. Remember that is the parameter A which should be compared, 
in rapidity intervals not too close to full phase space, with the average number of 
clans as obtained from NBD fits. One can clearly see in the figure the linear rise of 
A (Ay, W) as Ay grows, and the slow decrease with energy at fixed Ay. The result for 
nc{Ay, W) is reasonable but not fully satisfactory: in the limit Ay 0, nc{Ay, W) 
points to constant values, which differ at different energies, in contrast with the ex- 
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Fig. 3. Results of the GSPS model for the parameters of clan analysis as a function of the width 
of the rapidity interval j/c at different maximum allowed virtualities W = 50 GeV (solid line), 100 
GeV (dashed line) and 500 GeV (dotted line). Analytical solution with A = 2 and a = 1. 



pected ly-independent value ndO, W) = 1. The explanation for this anomaly lies in 
the approximation we had to use in the analytical calculations, which fail in the very 
small intervals Ay < 1. Indeed a Monte Carlo version of the model shows0 that the 
slope of X{Ay, W) when Ay matches that of n{Ay, W) in Fig. ||, so that actually 
one finds nc(0, W) = 1. 

6. Conclusions 

We have illustrated the GSPS model and its results; it is a parton shower model 
which was built from QCD-inspired splitting functions in virtuality and in rapidity, 
with Sudakov form factor normalization, to which the phenomenologically established 
idea of clans was added by allowing at each step in the cascading local violations of 
the energy-momentum conservation law (which is recovered globally in a statistical 
sense). The model has important predictive power in regions not presently accessible 
to full perturbative QCD; the results on clan analysis are consistent with what we 
know of single gluon jets disentangled using a jet finding algorithm from the Jetset 
Monte Carlo program and analyzed at parton level by assuming generalized local 
parton-hadron duality!. These predictions can (and hopefully will) be tested in the 
near future as clean samples of single gluon and quark jets of different energies have 
been separated at LEpi ; some preliminary data are available in full phase space only, 
an analysis of which is under wayS. 
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